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The nearby Supernova 1987A was accompanied by a burst of neutrino emission, which indicates 
that a compact object (a neutron star or black hole) was formed in the explosion. There has been 
no direct observation of this compact object. In this work, we observe the supernova remnant with 
JWST spectroscopy, finding narrow infrared emission lines of argon and sulfur. The line emission 

is spatially unresolved and blueshifted in velocity relative to the supernova rest frame. We interpret 
the lines as gas illuminated by a source of ionizing photons located close to the center of the 
expanding ejecta. Photoionization models show that the line ratios are consistent with ionization 

by a cooling neutron star or a pulsar wind nebula. The velocity shift could be evidence for a 


neutron star natal kick. 


upernova 1987A (SN 1987A) in the Large 
Magellanic Cloud is the most recent super- 
nova (SN) in the Local Group of galaxies, 
so it has been extensively studied. It origi- 
nated from the collapse of a blue supergiant 
star of 15 to 20 solar masses (Mo). Observations 
taken 35 years after the explosion show that the 
remnant of SN 1987A contains several distinct 
components (Fig. 1) (7). The inner ejecta, con- 
taining the synthesized heavy elements produced 
in the explosion, are now interacting with an 
equatorial ring (ER) and the surrounding cir- 
cumstellar medium (2). The ER is thought to 
have originated ~20,000 years before the ex- 
plosion, probably from a binary merger event 
experienced by the progenitor star (3, 4). The 
inner ejecta produce broad emission lines with 
velocities of 2000 to 4000 km s™, powered by 
the radioactive decay of “Ti nuclei and x-rays 
generated by the interaction of the ejecta with 
the ER (J). 
A compact object, such as a neutron star (NS) 
or black hole (BH), is expected to be located in 
the central region but has not been observed. A 


neutrino burst was detected coinciding with 
SN 1987A (5-7); its duration indicates that the 
compact object is probably a NS (8). Ionizing 
photons emitted by such a NS have been pre- 
dicted to excite emission lines from heavy 
elements in the ejecta, with velocities of 100 
to 1000 km s™ (9). In this study we searched 
for those emission lines. 


Infrared spectroscopy of SN 1987A 


We observed the SN 1987A remnant on 16 July 
2022 (day 12,927 after the SN explosion) using 
the JWST with the Mid-Infrared Instrument/ 
Medium Resolution Spectrometer (MIRI/MRS) 
and the Near-Infrared Spectrograph/Integral 
Field Unit (NIRSpec/IFU). Both instruments 
produce data cubes, which contain a spectrum 
for each spatial pixel (10). 

The data cubes exhibited an excess of emis- 
sion in the central pixels of the ejecta (Fig. 2) at 
observed peak wavelengths of 4.5293 + 0.0003 um 
in the NIRSpec data and at 6.98606 + 0.00003 um 
in the MRS data, with the latter line being much 
brighter. We identify these as forbidden emission 


4.52922 and 6.985274 um, respectively. Figure 2, A 
to P, shows velocity slices of the MRS data cube 
centered on the wavelength of the [Ar 1] line, and 
Fig. 2, Q and R, shows slices of the NIRSpec cube 
containing the [Ar vr] line. The [Ar 1] line in 
SN 1987A has been observed previously (11), 
but the spatial and spectral resolution were in- 
sufficient to separate the emission from the 
ejecta and the ER. The JWST observations show 
that it arises from a central source, separated 
from the ER. 

Figure 3A shows the [Ar 11] line profile from 
the central source, which has a narrow strong 
peak and a weaker broad extension to the red. 
We fitted this asymmetric line profile with two 
Gaussian functions. Figure S2 compares the 
line from the center of the ejecta with the cor- 
responding [Ar 1] line from the ER. The narrow 
emission line component in the center is dis- 
placed by -259.6 km s™ with respect to the rest 
frame of SN 1987A, which has a heliocentric 
velocity of +286.5 km s™ (12), whereas in the 
spectrum of the entire ER, the [Ar 1] line shows 
a much smaller blueshift of -12.2 + 1.6 kms + 
(10). The velocity difference and spatial dis- 
tribution indicate that the emission line from 
the central source is not scattered light from 
the ER. 

The MIRI spectral resolution at the [Ar 1] 
wavelength corresponds to a full width at half 
maximum (FWHM) of 87.7 + 1.1 km s™ (13). 
We measured that the narrower ejecta line com- 
ponent has a FWHM of 121.9 + 1.3 km s~}, so it 
is marginally spectrally resolved. The emission 
from the entire ER has a FWHM of 286 km s”, 
and the broader component has a value of 
362.4 km s™". Table 1 lists the velocity offsets, 
FWHMs, and luminosities that we measured 
for the emission lines from the central ejecta 
[using a distance to SN 1987A of 49.6 kpc 
a4)]. 

The NIRSpec spectral resolution at the [Ar v1] 
4.529-.m line is ~261 km s™. Fitting a similar 
two-component model as for the [Ar 1] line, we 
found that the peak of the [Ar vi] line has 
consistent velocity to that of [Ar 11], but the 
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broad components differ from the narrow 
blueshifted component (Table 1). Both the 
[Ar 1] and [Ar vr] peaks are consistent with 
being spatially unresolved [the spatial reso- 
lution FWHM is 0.350 and 0.212 for [Ar 11] 
and [Ar vi], respectively, where " indicates arc 
seconds (10)]. The positions of maximum emis- 
sion in the [Ar n] and [Ar vi] lines spatially 
coincide (within the uncertainties) with each 
other and with the previously determined cen- 
ter of the ER (Fig. 1), which is assumed to be the 
explosion site (15). 

In the MIRI data, we also detected lines of 
[S rv] at 10.51 um and [S m] at 18.71 um (fig. S3). 
Both sulfur lines show weak components from 
the central ejecta, displaced from the SN 1987A 
velocity by a similar amount as the argon lines. 
We can only place an upper limit on the lumi- 
nosity of an [Ar m] line at 8.991 um (Table 1). 
In the NIRSpec data, we detected the [Ca rv] 
line at 3.207 um (fig. S5); however, its velocity 
is shifted to the red and its maximum emission 
is located slightly north of the [Ar vi] peak, so it 
might be unrelated to the argon lines (70). 


Interpretation of the central emission 


The only narrow emission components from 
the ejecta that we identified are from Ar and 
S, which are both produced by the nuclear 
burning (fusion) of O and Si. This indicates 
that the emission arises from the inner core of 
the ejecta (fig. S7). Explosion models of SN 1987A 
(16) show that the innermost ~0.07Mo of ejecta 
is dominated by He and “Ni (the latter decays 
into Fe). The next zone outward, with a mass 
of ~0.3M.,, contains material produced by 
oxygen burning—Si, S, Ar, Ca, and, in its outer 
region, also O (10). 

The same elements have been observed in 
other young oxygen-rich supernova remnants 
(SNRs) from massive star explosions, including 
Cas A, SNR G54.1+0.3, and SNR 0540-69.3 
(17-20) (supplementary text). Both SNR G54.1+0.3 
and SNR 0540-69.3 contain pulsar wind nebu- 
lae (PWNe), powered by the spin-down power 
from a rapidly rotating NS with a strong mag- 
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netic field. By contrast, Cas A contains a central 
cooling NS (CNS) with a surface temperature 
~1.8 x 10° K (27) without a PWN. The emission 
lines in Cas A are mainly excited by a shock in 
the outer parts of the SN ejecta. 

The narrow [Ar 1], [Ar vi], [S n1], and [S 1v] 
lines that we observe from the center of SN 1987A 
must be excited by a source of ionizing photons 
or a shock wave. Potential sources include pho- 
tons from a PWN generated by a NS, photons 
directly from a CNS, photons emitted by accre- 
tion onto a compact object, or shock excitation 
by a PWN shock. We also considered five other 
possibilities (supplementary text): ionization by 
radioactive “Ti; excitation produced by x-rays 
from the ejecta and circumstellar medium in- 
teraction; a surviving companion star; re- 
flection of the narrow line emission from the 
ER collision by dust; or emission from an in- 
going, reverse shock. All five possibilities were 
excluded. 

The remaining potential explanations all in- 
clude the presence of a young NS or a BH in the 


Fig. 1. Optical image of 

SN 1987A taken in 2022, 

35 years after explosion. 
Hubble Space Telescope (HST) 
image in the F625W filter (10), 
which is dominated by Ha 
emission. The freely expanding 
inner ejecta and ER are 
labeled. White contours mark 
the [Ar vi] 4.529-um emission 
observed with NIRSpec (at 

40 and 70% of the maximum 
surface brightness in Fig. 2S). 
The white star denotes the 
center of the ER (15). 


Count rate (counts s-1) 


center. We consider a BH unlikely because the 
progenitor star of SN 1987A is thought to have 
had a total mass of < 20M. and an Fe core mass 
of $ 2Mọ (22). The explosion ejected at least 
0.07Mo of 56Ni (23). That leaves less core 
mass remaining than would be required to form 
a BH, which is = 2.2Mọ (24). Accretion of ejecta 
bound to the NS could be a possibility. The 
estimated luminosity from this, ~2 x 10” erg s” 
at 35 years (25), is too low compared with the 
observed luminosity (Table 1), but the relevant 
physical parameters are uncertain. 

A hot CNS would have a temperature of 
> 10°K and a luminosity of = 10°*ergs~! (26), 
which is sufficient to ionize part of the inner 
ejecta (10). An analogous case is the NS in Cas 
A, although it is older and fainter than what 
we expect for SN 1987A. 

The energy loss from a pulsar is in the form of 
electrons, positrons, and magnetic fields, which 
produce an expanding high-pressure PWN. 
High-frequency synchrotron radiation from the 
relativistic particles can then ionize the ejecta 


Table 1. Parameters for emission line components from the central ejecta. Luminosity calculations adopt a distance to SN 1987A of 49.6 kpc (14). No 
broad component was required for the [Ar m], [S v], and [S i] lines (fig. S3). The [Ar vi] lines were measured in the NIRSpec data; all other lines are from the 
MIRI data. Uncertainties are lo errors on the fitting parameters; upper limits are 30. Dashes indicate not applicable. 


Narrow Broad 
component component 
Line Wavelength (um) 
Velocity FWHM Luminosity Velocity FWHM Luminosity 
offset (km s~’) (km s~) (10°° erg s’) offset (km s~) (km s~) (10°° erg s’) 
LAO eee ee a ae eee SS eee, Ie oe eee ee Seve ae 2.88 + 0.40 | 
[Ar u] 6.985 -259.6 + 0.4 WAL) ge ils) I'S) ae DU -200.6 + 3.7 362.4 + 8.0 104.1 + 3.0 
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Fig. 2. JWST observations of the [Ar 1] and [Ar vi] lines in SN 1987A. 

(A to O) Velocity slices of the MIRI/MRS data cube around the [Ar 1] line at 
6.985 um. Labels in each panel indicate the corresponding velocity ranges, 
measured with respect to the systemic velocity of SN 1987A (12). Emission is 
apparent from a central source, which has its maximum intensity offset in 
velocity from the ER by ~-250 km s™. The white circle in (F) indicates the region 


(9). The luminosities and spectra of the PWNs 
discussed above are not consistent with the 
bolometric and monochromatic luminosity 
of the SN 1987A remnant, whereas a slowly 
rotating and/or weakly magnetized NS is 


compatible with the observations (15). We 
expect an expanding PWN bubble to sweep 
up adense shell of ejecta, separated from the 
PWN by a shock, which could also excite 
the lines that we observed. The velocity of the 
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D [Aru] -444 > -377 km s~ 


H [Ari] -170 — -101 km s- 


Surface brightness (MJy sr-1) 


Surface brightness (MJy sr-") 


used to extract the spectrum in Fig. 3A. (P) Residuals between (F) and (A). 
(Q) NIRSpec observations of the [Ar vi] 4.529-um line, combining all velocity 
slices between -450 and +380 km s 1. (R) The NIRSpec data stacked outside that 
velocity range. (S) Residuals between (Q) and (R), separating the emission from 
the central source. The bright point in the lower left areas of (Q) and (R) is an 
unrelated star. Color bars are in units of megajansky (MJy) per steradian. 


dense shell would depend on the energy input 
from the pulsar along with the density and 
expansion velocity of the ejecta (9). These 
parameters are all uncertain, but estimates 
for SN 1987A produce a shell velocity of 150 
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Fig. 3. Velocity profiles of the [Ar 1] 6.985-um and [Ar vi] 4.529-ym 
lines from the central source. (A) The [Ar 11] line in the MRS data (black 


histogram) extracted from the region indicated in 


is a model fitted to the data, consisting of two Gaussian components 


to 300 kms‘ and a PWN shock velocity of 30 
to 50 km s™ (10). 


Photoionization models 


To investigate these interpretations in more 
detail, we updated a photoionization modeling 
code (9) and applied it to the PWN and CNS 
scenarios (10). The resulting emission line spec- 
trum is primarily determined by the ionization 
parameter © = Lion /Mion?”, Where Lion is the 
ionizing luminosity, Mion is the number density 
of ions, and r is the distance to the ionizing 
source; it is affected to a lesser extent by the 
spectral form of the luminosity L(v), where v is 
the frequency (10). For the ionizing spectrum 
of a PWN, we assume a power law L(v) © v °, 
with o = 1.1. This matches the ultraviolet-to- 
x-ray spectrum of the Crab Nebula (27) and is 
within the range of 0.5 < a < 1.1 previously 
found for x-rays in young PWNs (28, 29). We 
also investigated using steeper power laws 
and found that it has little effect (supple- 
mentary text). Previous work has constrained 
Lion S 10% erg s~ at energies above 13.6 eV (15). 

For the young CNS scenario, previous models 
have found an upper limit to the thermal lu- 
minosity of <4 x 10% erg s~! and a surface tem- 
perature of ~(1.5 to 3) x 10° K at 35 years after 
the explosion (26, 30, 31). Central compact ob- 
jects, like the NS in Cas A, are thought to have 
carbon-dominated atmospheres (21), although 
alternative models producing a hard spec- 
trum have been discussed (32). We adopted a 
carbon-dominated CNS spectrum (27) and set 
Lion = 3 x 10” ergs7t. 

The models depend on the density in the 
ejecta core, which is uncertain. After explosion, 
the mass inside the oxygen core of a 19—Mo5 
progenitor is ~3.0M., after excluding the 
material that forms the NS (16). The line pro- 
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Fig. 2F. The orange curve 


files indicate a core velocity of ~2200 km s™ 
(0). For a uniform core density, Mion 7 2.6 x 107 
(4/22) cm”? after 35 years, where A is the 
mean atomic weight in the O-Si-S-Ar core. 
However, both observations (33-36) and simu- 
lations (37) indicate that clumping and fila- 
ments produce regions with higher densities 
in this zone, whereas Fe-rich regions can have 
lower density as a result of early heating by 
decay of °°Ni (10). If a PWN is present, it pro- 
duces a dense shell (as discussed above), which 
is likely to be unstable and clumpy (9, 38). 
We therefore assume a density of Nion = 2.6 x 
10* cm~’, which is 10 times the average esti- 
mated above, in both models. 

Figure 4 shows the lines predicted by the 
models of the PWN and CNS scenarios for an 
ejecta composition (fig. S7) typical of explosive 
oxygen burning in a SN with a 19—Mo proge- 
nitor star (10, 16). The blueshifted line veloci- 
ties and dust effects indicate that only a fraction 
of the total ionizing luminosity is converted 
to the observed emission lines (supplemen- 
tary text). We therefore normalized the [Ar 1] 
luminosity by a covering factor, chosen to 
match the observed luminosity, and scaled 
the other lines to their relative model lumi- 
nosities (70). For consistency, we included 
only the luminosities of the blueshifted nar- 
row component. 

The relative luminosities of the [Ar 1] 6.985-um 
and [Ar v1] 4.529-um lines are reproduced in both 
models (Fig. 4), with this assumed density and 
E = 0.3 (10). Both models predicted other high- 
ionization emission lines, including [Ca 1v] 
3.20741m, [Ca v] 4159-um, [Ar v] 7.902-um, [Ar m] 
8.991-um, [S 1v] 10.51-um, [Ar v] 13.10-um, [S m] 
18.71-um, [Ar m] 21.83-um, [O 1v] 25.89-um, and 
[Fe n] 25.99-um lines (Fig. 4). Most of these 
lines have previously been observed in young, 


23 February 2024 


-500 0 
Velocity (km s71) 


500 1000 


(blue and gray dashed curves). (B) The same [Ar 11] data (black histogram), 
compared with the [Ar vi] line from the NIRSpec data (red histogram). 

The [Ar vi] data have a lower spectral resolution and have been scaled by a 
factor of 76 for display. 


oxygen-rich SNRs (supplementary text). How- 
ever, in SN 1987A, the lines at wavelengths higher 
than ~8 um are much weaker than predicted 
(relative to the [Ar 1] line) or are absent en- 
tirely, although we do detect faint emission for 
[S rv] 10.51-um and [S m] 18.71-um lines (fig. S3). 


Dust effects on the observed emission 


Infrared observations have shown that a large 
amount of dust was formed in the ejecta of 
SN 1987A (39, 40), consisting of a mixture of 
grains composed of silicates and graphite (41, 42). 
However, the composition of the dust close to the 
center is likely to differ from that of the outer 
regions (10). Dust nucleation models have pre- 
dicted that silicates dominate over carbon-rich 
dust in the central regions (43-46), of which 
Mg,SiO, (forsterite) is expected to have the 
highest abundance. Observations of young SNRs 
have shown that silicates are the dominant type 
of dust (0). 

Silicate dust grains strongly absorb wavelengths 
longer than ~8 um (fig. S8). This is a potential 
explanation for why we do not detect the long- 
wavelength lines at the luminosities predicted 
by the photoionization models. 

We are unable to estimate the optical depth 
Tabs at those wavelengths because of the un- 
known dust composition, grain size, and clump- 
ing. However, Taps is likely to be large (47). Figure 
4 also shows the dust absorption expected for an 
assumed optical depth of Tabs = 6.5 at 10 um and 
amorphous forsterite composition (48) and 
the line luminosities after including this dust 
absorption. Although the predicted luminosi- 
ties of lines shorter than ~8 um are only slightly 
affected, those at longer wavelengths are sub- 
stantially suppressed. For [S 1v] 10.51-um and 
[S mr] 18.71-um lines, this correction brings 
the model predictions close to the observed line 
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Fig. 4. Photoionization models compared with observed lines. (A) Model of the 
oxygen-burning layers expected in a SN ionized by a young CNS. (B) Model of 
the same material ionized by a PWN. The optical depth due to silicate dust 
(magenta curve, on the right axis) assumes Mg2SiO4 composition and taps = 

6.5 at 10 um in both cases. Circles show the model-predicted line luminosities (left 


luminosities. The observed upper limit on 
the [Ar n] 8.991-um line is consistent with 
the model. 

The CNS model (Fig. 4A) predicts strong 
[Si 1] lines at 1.607 and 1.645 um. These are not 
detected in our observations, perhaps because 
they are blended with broad lines from the 
ejecta. Dust scattering, dust depletion, and ioni- 
zation of silicon by the ultraviolet flux from the 
ejecta and ring could all reduce the predicted 
luminosities of these lines (170). 

We have also modeled the expected line 
luminosities from a PWN shock (10). Most of the 
observed lines are reproduced (supplementary 
text) (fig. S12) without requiring dust absorp- 
tion. However, the [S m] 18.71-um line is under- 
produced by a factor of ~25. 


The nature of the compact object 


Given the uncertain model assumptions about 
the ionizing spectrum, density, elemental abun- 
dances, and dust properties, we conclude that 
the observed lines can be explained by either 
the CNS or PWN photoionization models. Both 
require the presence of a NS. Because thermal 
emission from the CNS should always be pres- 
ent at some level, a combination of both models is 
also possible. For the model parameters adopted 
above, the PWN shock scenario is less likely but 
not excluded (70). 

There have been previous indications of 
emission from a central object in SN 1987A. 
Submillimeter dust maps of the ejecta show a 
peak in the temperature close to the central 
region, which might be due to extra energy 
input from a central source (47). The peak tem- 
perature position does not coincide with the 
location of maximum [Ar vi] emission in our 
observations (10); however, this does not rule 
out additional heating from a compact object 
(31). X-ray observations of SN 1987A have de- 
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than the symbol size. The 


tected nonthermal emission at energies above 
10 keV, which has been interpreted as being 
due to a central source (49, 50); this would 
support the PWN scenario. However, the same 
observations have alternatively been interpreted 
as thermal emission from the shocks in the 
ER (51). 

The position and velocity of the argon emis- 
sion lines constrain any natal kick—a velocity 
imparted on the NS during its formation by 
anisotropies in the explosion mechanism. The 
offset in position of the narrow emission com- 
ponents corresponds to a velocity 416 + 206 km s* 
given the time since the explosion (supplementary 
text). However, converting this to a kick velocity 
depends on the separation between the NS and 
the line-emitting region, which differs between 
the CNS and PWN models. 

Our observations of high-ionization emission 
lines in the center of the SN 1987A remnant 
indicate the presence of a source of ionizing 
photons, probably a NS. We cannot determine 
which scenario is more likely—young CNS or 
PWN-—but both require the presence of a NS. 
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